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ABSTRACT 

In this research work, the Klein-Gordon equation for Modified Eckart plus Inverse 
square potentials using the generalized parametric form of Nikiforov-Uvarov method 
was studied. Energy eigenvalues and the corresponding unnormalized wave function 
expressed in terms of Jacobi polynomial was obtained. It was discovered that the 
two special cases of this potential comprises of the Eckart potential and the Inverse 
Square potential. 

Keywords: Parametric Nikiforov-Uvarov method, new approximation scheme, 
Eckart, plus Inverse square potential. 


1.0 INTRODUCTION 

In the last decade, a lot of articles have focused on the analytical and 
approximate solutions of the Schrodinger, Klein-Gordon, and Dirac equations. This 
is because the solutions of these wave equations contain all the necessary 
information for the quantum system under consideration. As is well known, the 
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Schrodinger wave equation is used to describe non-relativistic spinless particles, the 
Klein-Gordon equation for relativistic spin-zero particles, the Duffin-Kemmer-Petiau 
equation for relativistic spin-zero and spin-one particles, and the Dirac equation for 
relativistic spin-half particles. Many potential models have been studied using these 
wave equations, such as the Kratzer potential, Morse potential, Rosen-Morse 
potential, Manning-Rosen potential, Deng-Fan potential, and Poschl-Teller potential 
among others. See for equations [1-8] for details. This paper will discuss the 
relativistic effects, when a particle moves in the potential consisted Modified Echart 
plus Inverse square potential (MEISP) which is given as: 


Vis) = 


VqS a Co ac ± ac 2 

(l-S) 2 K 2 K 2 (l-S) K 2 (l-S) 2 


( 1 ) 


Solutions of bound states of the Klein-Gordon equation with vector potentials are 
obtained, the normalized angular wave function and radial wave function are 
presented, simultaneously, energy spectrum equations are also obtained. The paper 
is organized as follows: Section 1 has the introduction, the NU method is reviewed 
in section 2. In section 3, the Klein-Gordon equation with MEISP is solved using the 
NU method. Finally, we give a brief discussion in section 4 before the conclusion in 
section 5 and then the references also given. 


2.0 GENERALIZATION OF NIKIFOROV-UVAROV METHOD 

The main equation which is closely associated with the method is given in the 
following form (Nikiforov and Uvarov, 1988). 

r( S ) + ^ns)+j%iP(s) = o ( 2 ) 

Where a (s) and a(s) are polynomials at most second-degree, T(s) is a first-degree 
polynomial and <//(s) is a function of the hypergeometric-type. 

In order to find the exact solution to Eq. (2), we set the wave function as: 

ip(x) = 0(s)X (s) (3) 

And on substituting Eq. (3) into Eq. (2), then Eq. (3) reduces to hypergeometric- 
type, 

a(S)X"(s) + t (S)X'(S) + AX (5) = 0 (4) 

Where the wave function 0 (s) is defined as the logarithmic derivative 
0'OO = 7 r(s) 

000 <r(s)' 1 J 
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Where n (s) is at most first-order polynomial. 

The hypergeometric-type function 0 (s) polynomial solutions are given by the 
Rodrigues relation 

^ (s) = ^i [CT " (s)p(s)] (6) 

Where B n is the Normalization constant and the weight function p(s) most satisfy 
the condition 

^[o- n (s)p(s)] = T(s)p(s) (7) 

Where; 

t(s) = T(s) -I- 2n(s) (8) 


In order to accomplish the condition imposed on the weight function p(s) , it is 
necessary that the classical or polynomials f( s) be equal to zero to some point of 
an interval (a , b ) and its derivative at this interval at o ( s)>0 will be negative, that 
is; 

dr(s) 
c[s < 0 


Therefore, the function 7r(s) and the parameter A required for the NU method are 
defined as follows: 


tt(s) = 


A 


'O' — x> 


— o + Ko 


(10) 


Where A=k+n' (s) 


The parametric generalization of the NU method is given by the generalized 
hypergeometric-type equation as 


roo + 


~ (Ci-C 2 S) 

s(i - c 3 s ) 


ip'(S) + 


-^S 2 + f 2 5 - 

s^i-c 3 sy 


ip^s) = 0 


(11) 


Equation (11) is solved by comparing it with Eq. (2) and the following polynomials 
are obtained: 


T(i) = (Ci - C 2 S),ff(s) = 5(1 - C 3 S),a(s) = s 2 + - f 3 (12) 

Now substituting Eq. (12) into Eq. (11), we find 
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<r(s)C 4 + C 5 S + J[(C 6 - C 3 tf ± )S 2 + (C 7 + K ± )S + C 8 ] 


( 13 ) 


Where 

c 4 = \ (1 - c t ), C 5 = i (C 2 - 2 C 3 ), C 6 = C| + fj, C 7 = 2 C 4 C S - f 2 , C 8 = C 4 2 + 
^3, Cg = C 3 C 7 + C 3 C 8 + C 6 C 10 = C 4 + 2C 4 + 2V” Cg C 4 (14) 

The resulting value of k in Eq. (13) is obtained from the condition that the function 
under the square root be square of a polynomials and it yields, 


K± — ( C 7 + 2 C 3 C 8 ) + 2 yj Ccj C Q 

Where C 9 = C 3 C 7 + C 3 C 8 + C 6 

The new n (s) for k becomes 

7r(s) = C 4 + C 8 S — [(n/Cg + C 3y [C~ 8 )S — yfCo\ 

Using Eq. (8), we obtain 

t(s) = C 1 + 2C 4 — (C 2 ~ 2 C 5 )s — 2 [(-^Cg + C 3 yj~C^)s — 


(15) 


(16) 


(17) 


We obtain the energy equation as 

(C 2 — C 3 )n + C 3 n 2 — (2 n + 1)C 5 + (2 n + 1 )( A /C^ + C 3 J~C^) + C 7 + 2C 3 C 8 + 2yj C 8 C 9 = 0 (18) 


While the wave function is given as 


c 13 


u 13 / Cll \ 

Ms) = N ni s c ^( 1 - c 3 s) C3 P n - 2 C 3 S) 


(19) 


Where P n is the orthogonal polynomials. 

Given that = P n <“■« = E?-o , r ( n + a+1 Mn+g+i) fci) r (£±1) n-r ( 2 o) 

n ijr - l) r(fl+r+l)r(n+p-r+l)(n-r)!r! V 2 ) V 2 ) v J 

This can also be expressed in terms of the Rodriguez's formula 

1 . d \ n 

P n M \x) = —{x- !)““(* + lW—J (lx - l) n+a (x + 1)" + P) (21) 


3 SOLUTIONS OF THE RADIAL KLEIN-GORDON EQUATION 

The Radial Klein-Gordon Equation with vector V(r), potential in atomic units (h = c 
= 1) is given as; 
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d 2 R(r ) 
dr 2 


+ [{E 2 - M 2 ) - (E + M)V(r)\ R(r ) = 0 


( 22 ) 


The Modified Eckart Potential is given as: 

V(r) = ( V ° e ~ ar ) 
v J \(i -e- a n 2 J 


(23) 


The Inverse Square Potential: V(r) = — 


(24) 


Applying the transformation S — e ar , the new improved approximation is given as 

,-ar \ 


1 _ 1 
r 2 k 2 


^(^) + c 2 (^y 


(25) 


( i 4.0 — cir b \ “ 

^r) >A = 2(* ar = iy 


Where 


The superposed potential can be represented as MEISP 

j/(s) — v ° s _|_ AC ° _|_ ACl _|_ AC2 


(i-s ) 2 k 2 a: 2 (i-s) K 2 (i-sp 


(26) 


cL ^ i? ( 5 ) 

Applying the new improved approximation and after lengthy algebra, we have + 

^ s d T + (djv [W 2 - H + 2 + +Q-P + ms + (2ft 2 - H - Q - !]R(S) =0 


(27) 


Where 


2 = (e*+mI\ = f_ji\ AC 

r V 4a 2 J v \2 k 2 a 2 h 2 J 1 


H 


= (v^) Ac o. p = (iK/ = feraK. 






< 2 a 2 h 2 ' 


<2 a 2 h 2 


2 k 2 a 2 h 2 ' 


-1 1 

Ci = c 2 = c 3 = 1, c 4 = 0, c 5 = — , c 6 = - + 2/? 2 - H, c 7 = -4/? 2 + Q - P - 2 H, 


c 8 — 2/? 2 H Q J,c 9 - P J, c 10 — 1 + 2 


N 


2B 2 - H -Q - J,c 11 = 2 + 2\ 


-~P-] + ^2f 2 -H-Q-J 

(28) 
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Using the eigenvalue equation, the energy eigen spectrum of MEISP is found to be 

2 


P 2 = 


(P + Q + 2J) — (n 2 + n — ^ 

) - (2 n + 1)^| 

4 r 7 

(2 n + 1) + 2 J 



-( H + Q+J ) 


(29) 


The above equation can be solved explicitly and the energy eigen spectrum of MEISP 
becomes 


E 2 - 


M 2 —4a 2 




(30) 


4. EIGEN FUNCTION CONSIDERATION 

The wave function p(s) is given as 

^■11 * 

p(s ) = S^o^Cl - C 3 S)~ Cl0 ~ 

Using equation (28), we get the weight function as 

p( S ) = s u (i-sr 

Where U = p 2 -R-P-X, and V = 2 J^~P~ J 

Also we obtain the wave function X n (s ) as 
X„(s) = P^\l - 2s) 

Where pf’ v ' are Jacobi polynomials 

Lastly, cp(s) = s Cl2 ( 1 - C 3 s)~ Cl2 
Using equation (28) we get 

u 1+v 

cp(s) = 52(1 — 2)~ 


(31) 


(32) 


(33) 


(34) 


(35) 


We obtain Radial wave function from equation (28) 

R n (s) = N n cp(s)X n (s ) (36) 
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As 


u l+v ( ^ 

«n(s) = JV„Sl(l - s)— P n (U - V) ( 1 - 2s) (37) 

Where n is a positive integer and N n iS the normalization constant. 

5. DISCUSSION 

In this subsection, we will consider some special cases of the potential in 
consideration. 

(i) If we set V 0 - 0 in equation (23), the eigen energy spectrum of equation 
reduces to (34) Which is the energy spectrum for inverse square potential 
using the new approximation scheme. 


E z — M 2 — -4a 


2^ 


jv2 K^a^hZy 




(34) 


(2n + l)+J-(^)AC 2 


J J 


(ii) If we set A - 0 in equation (24), the eigen energy spectrum of equation 
(29) reduces to (35) Which is the energy spectrum for modified Echart 
potential under the new approximation scheme 


E 2 -M 2 = -4a 


i(l^) v <>)< n2+n -l)^ 2n+1 ^^) v < >] 

(2n + l 


(35) 


6. NUMERICAL ANALYSIS 

In this subsection, we discussed the behavior of our calculated Eigen energy values 
and the superposed potential (MEISP) at different values of the screening 
parameters. Figure I shows the variation of the eigen energy with screening 
parameter. It can be observed that eigen energy increases with increasing screening 
parameter. Figure II shows the relationship between the superposed potential with 
the internuclear distance. The plot has a repulsive behavior within the potential 
energy diagram. Table I and II depicts the numerical result of the eigen energy 
and superposed potential at different screening parameters. 
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Variation of Energy with Principal Quantum 
Number 



n 


Figure I: Plot of eigen energy versus Principal Quantum Number 


MEIS Potential versus Internuclear distance 



Figure II: Variation of Potential with internuclear distance 


Table 6.1: Eigen energy at different Screening Parameter 


E, 

n 

o = 0. 01 

o = 0. 02 

1 

15.64899551 

15.71069848 

2 

15.68466633 

15.85621798 

3 

15.73857948 

16.08133291 
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4 

15.81123382 

16.39499224 

5 

15.90331266 

16.81052189 


Table 6.2: Potential Values at different Screening Parameter 


R(um) 

V(r) 


o = 0. 01 

o = 0. 02 

0.01 

105020002.30 

262660002.26 

0.02 

26260002.260 

6567502.2630 

0.03 

11673335.600 

2920002.2400 

0.04 

6567502.2630 

1643126.8520 

0.05 

4204002.2630 

1052001.8420 


7. CONCLUSION 

In this paper, we investigated the solutions of the Klein-Gordon equation for the 
Modified Eckart plus Inverse square potentials using a new approximation method 
for the centrifugal term via parametric Nikiforv-Uvarov method. The eigenvalues of 
the potential reduced to that of well-known potentials viz Echart potential in (35) 
and Inverse square potential in (34), when we make appropriate choices of the 
parameters. Finally, we also obtained the numerical result of our eigen energy 
spectrum. 
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